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Abstract 

After reviewing the Hermitian one matrix model, we will give a brief 
introduction to the Hermitian two matrix model and present a summary 
of some recent results on the asymptotic behavior of the two matrix model 
with a quartic potential. In particular, we will discuss a limiting kernel 
in the quartic/quadratic case that is constructed out of a 4 x 4 Riemann- 
Hilbert problem related to Painlevé H equation. Also an open problem 
will be presented. 



1 Introduction 

The study of the local statistics in larga random systems of interacting partidas, 
such as the aiganvalues of large randoni matrices, is a central thenie in random 
matrix theory. The universality principie states that the local statistics obey 
laws that have a universal character and do not depend on the precise definition 
of the underlying probability measure but rather on some general characteristics. 
Generic examples are the sine universality in the bulk and Airy universality near 
a soft edge. However, there may be singular points near which the limiting local 
correlations are described by more complicated universality classes. In some of 
these cases remarkable connections to Painlevé equations have been found. 

A model on which substantial progress has been made is the Hermitian one 
matrix model (also called the Unitary Ensembles) . A milestone was the rigorous 
verification of the bulk and soft edge universality conjectures in P^[2Ü] based on 
the orthogonal polynomial approach and Riemann-Hilbert techniques. Among 
further important developments are the treatment of the three types of singular 
points that can occur. Near two of those the local correlations are described by 
Vt-functions associated to special solutions to members of the Painlevé I and II 
hicrarchies [3 [H Hü [H [S5] . 

In this paper, our main focus will be on the Hermitian two matrix model 
which is a natural extensión of the Hermitian one matrix model. There is strong 
evidence that the two matrix model gives rise to a family of singular situations 
that is much larger than the one for the one matrix model. For example, it is 
believed to be a good model for generating the (p, q) conformal minimal models 
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[T7] . The classification and characterization of the local correlations near the 
singular points is an important open problem in random matrix theory. 

By using the Riemann-Hilbert approach for the associated biorthogonal poly- 
nomials, we recently analyzed the asymptotic behavior of the two matrix model 
with one quartic potential [331 USl [23 HHI- Here we will report on that 
progress. As it is important to compare our results with the results for the one 
matrix model, we provide in Section 2 a brief overview of the one matrix model 
and the singular cases that can appear. In particular, we review the results on 
the case of one to two interval transition related to the Painlevé II equation. 
This situation will also return in our discussion on the two matrix model. 

In Section 3 we discuss the definition of the Hermitian two matrix model in its 
general form and the relation with certain biorthogonal polynomials. In Section 
4 we first discuss the vector equilibrium problem that was a key ingredient 
for the asymptotic analysis [351 123 HH] of the Riemann-Hilbert problem for 
the biorthogonal polynomials. Second, we present a phase diagram and a new 
critical phenomenon for the quartic/quadratic case that we analyzed in ¡22] . 
Interestingly, the limiting kernel is constructed out of a 4 x 4 Riemann-Hilbert 
problem associated to the Hastings-McLeod solution to the Painlevé II equation. 

Finally, in Section 5 we show that the singular case of [33] can be embedded 
into a larger class of singular cases for which a rigorous analysis may be within 
reach. We leave this as an open problem. 

2 One matrix model 

The Hermitian one matrix model is defined as the probability measure on the 
space of n X n Hermitian matrices given by 

(2.1) ^exp(-nTrl/(M))dM, 

where V is a polynomial of even degree and positive leading coefñcient, dM = 
nr=i "i-^" ni<j 'if'-^-^ii dImMy, and Z„ is a normalizing constant. We will 
be concerned with the limiting behavior of the eigenvalues of a matrix M taken 
randomly from (|2.ip as n — >■ oo. 

2.1 Global limit 

By integrating over the unitary group the probability measure (|2.ip reduces to 
the following probability measure on the eigenvalues 

(2.2) -i- ]^(x, - a;j)2e-"^^=i ^^^^^ dxi • • - dx™. 

Because of the Vandermonde determinant the probability of for eigenvalues 
to cluster is small and henee the eigenvalues appear to repel each other. On 
the other hand, due to the exponential factor the probability of finding the 
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eigenvalues spread out is small and henee V acts a confining potential. As 
n — >■ oo the eigenvalues seek the equilibrium situation for these two competing 
mechanisms. More precisely, the empirical eigenvalue distribution i 
converges weakly to almost surely, where is the unique minimizer of the 
energy functional 

(2.3) iviiy) = JJ logJ-^^d^l{x)d^i{y) + J v{x)d^l{x), 

minimized over all probabihty measure /i on M. 

This minimization problem can be solved in terms of an algébrale curve. 
There exists a polynomial Q oí degree deg V — 2 such that the function 

= V'{x) - I -^áMs), 



is a solution to the equation 

(2.4) e - V'{x)Í + Q{x) = 0. 

This curve can be derived using the variational conditions corresponding to 
minimizing (j2.3p . It can also be derived by a change of variables in (I2.ip (see 
for example ¿39j)- It particularly implies that /iy is absolutely continuous with 
a density 

d/iy , , 1 



ax IT 

where is the negative part of R{x) ~ V {x)'^ /A— Q{x) (so R = R+ — R-). In 
particular, one readily verifies that the support S{iJ,v) of /iy consists of a finite 
number of intervals. In the generic situation, the density is strictly positive in 
the interior and vanishes as a square root at the endpoints. 



2.2 Airy and sine universality 

Although the eigenvalue distribution on the global scale has a deterministic 
limit, interesting point processes are observed by scaling the eigenvalues around 
a point in the support of fiy such that the average distance between eigenvalues 
is of finite order. Before we describe the asymptotic behavior on the local scale, 
we discuss the orthogonal polynomials that intégrate the one matrix model. 

For k £ Nq let pk^n be the unique monic orthogonal polynomial of degree k 
with respect to e~"^(^)da; on R, i.e. 

y"pfc,„(x)x^e-"^("Mx = 0, j = O, . . . , fc - 1. 

Moreover, let Kn be the reproducing kernel 

n — 1 

K^{x,y) = e-"(^(-)+^(^))/2 J2 T^pU-)Pk,n{yl 
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where h\^ — \\Pk,n\\1^,^(^^-r^v{:,)¿^y Then the eigenvalues of M taken randomly 
with respect to (|2.1I) form a determinantal point process with kernel -fí„. This 
means that (|2.2p and the marginal densities (or, up to a scalar, the correlation 
functions) can be written as 

(2.5) j ' " j ^(^1' ■ • -i^yMxk+i ■ ■ -áxn = det {Kn{xi,Xj))\ -^-^ , 

n — k times 

for k = 1, . . . ,n. For more details on determinantal point processes we refer 
to [TÜl ISHl Sni US EZl ÜS- The main point is that íí„ characterizes the point 
process entirely and in order to find the asymptotic behavior of the process on 
the local scale it sufRces to compute the scaling limits of the kernel ií„. 

Let a:* G S{nv) such that p{x*) > 0. Then for large n the distance between 
eigenvalues near x* is of order ^ and we have the limit 

1 / ^ X ^ V \ sinTríx — v) 
(2.Q) lim K X H X + I — 

n^oc np{x*) "\ np{x*)^ np(x*) J Tr{x — y) 

uniformly for x,y m compact subsets. Note that the right-hand side is indepen- 
dent of X* and V. This is the celebrated principie of sine universality in random 
matrix theory. In [TÜ! HD] it was proved by using RH methods for analytic po- 
tentials V and later extended in [JSl [SI] to more general situations. It should 
also be noted that in recent years there has been an interesting development 
on proving (|2.6p by using more classical methods for orthogonal polynomials, 
obtaining even more general results. See ¡46] for a survey. 

Now let X* be an endpoint of an interval in the support of pv and assume 
that we are in the generic situation so that the density vanishes as a square root 
at X*. Then the distance between eigenvalues near x* is of order ~ and 
we have the scaling limit 
(2.7) 

r ±1 ^ . ^ ^ *, y \ Ai{x)Ai'{y)~Ai'ix)Ai{y) 
lim - — TTnTrK,, x ± - — TTn7T,x ± — 



(cn)^/-^ "\ (cn)^/-^ ' {cn)"^^^ J x — y 

with + is if X* is a right endpoint and — if x* is a left endpoint. Here Ai stands 
for the Airy function. 



2.3 Singular points 

The scaling limits (|2.6p and (|2.7p are valid for the regular points of the spectral 
curve (j2.4p . It may happen that there are singular points and around these 
points we obtain different and more complicated limits. There are three types 
of singular points that can occur. 

1. Exterior singular point a point outside the support where we have 
equality in the variational inequality for the equilibrium problem ([2.3^ : 
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2. Interior singular point a point x* in the interior of such that 
p{x) = c(x — a;*)^"'(l + 0(1)) a,s X ^ X* for some m G N and c > 0; 

3. Singular edge point an endpoint x* oí one of the intervals in ¿"^j, such 
that p{x) = c{x ± a;*)^™+-^/^(l + o(l)) as a; — )■ x* for some m e N and 
c > O (here ± depends whether we have a right or left endpoint). 

By varying the coefficients of V these singular points correspond to transitions 
in the support S{fj,v)- For example, we can let an interval in the support shrink 
to a point after which it disappears, leading to an exterior singular point. By 
letting two intervals merge to one we créate an interior singular point. If both 
transitions happen simultaneously near the same point we obtain a singular 
edge point. By involving several intervals simultaneously we obtain the higher 
order singularities. 

Near singular points, the local correlations have limits that are different from 
(|2.6I) and (I2.7P . The singular edge points and interior singular point are related 
to special Solutions to members of the Painlevé I and Painlevé II hierarchy 
respectively, which we will discuss in more detall. The treatment of exterior 
singular points of the first type can be found in [51 [T31 [SD] but they are not 
related to the Painlevé transcendents. 



2.3.1 Interior singular points 

We will now discuss the singular case of quadratic vanishing at an interior point 
in more detall. This situation will also play a role in the upcoming discussion 
on the two matrix model. 

Let US consider the double well potential 

V{x) = \x^-x^. 

In this case, the origin is a singular interior point and the density vanishes 
quadratically 

p(x) — ^x^\/2 — x^dx 

This situation was analyzed in [7 . For the treatment of more general potentials 
leading to quadratic vanishing at an interior point see [151 155| . 

The singular case is analyzed by means of a double scaling limit. We in- 
troduce a parameter A and consider the random matrix model with Vx — XV . 
For A < 1 the origin is a regular point and the 5*^^,- consists of one interval. 
For A > 1 a gap opens up at the origin and the support consists two intervals. 
For the critical valué A = 1 the density of the equilibrium measure vanishes 
quadratically at the origin. We take the limit n — > 00 but at the same time we 
let the gap open or cióse. That is, we let n — > 00 and A„ — !■ 1 simultaneously 
such that 

71^/3 (A„ - 1) ^ s e M. 

Then for some constants ci , C9 > O we have 
(2.8) 

1 / X y \ 
li'? , 7 üTí^" ^* + 7 ÜTT'^* + 7 ÜTT'^" ^ Kpu{x,y;c2s). 
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Figure 1: The jump contour in the complex C-plane and the constant jump 
matrices Jk on each of the rays Tjt, fc = 1, . . . , 4. 



Here Kpii{-, •; v) with e M is a family of kerncls for which the description is 
more compUcated than the limits in (|2.6p and (|2.7p . It can be characterized 
in terms of a Lax pair for the Hastings-McLeod solution to the Painlevé II 
equation. Here we wiU treat the Riemann-Hilbert problem(= RH problem) 
characterization of the kernel. See [33] for more details on the Riemann-Hilbert 
approach to the Painlevé equations. 

First we define the contour Yj<^ = Fi U r2 U U consisting of the rays 

Fi = e'''/^M+, Fa = e^'''/*^R+, F3 = -Fi, F4 = -F2. 

All rays are oriented towards infinity. The orientation also allows us to provide 
each point of the contour with a + side lying at the Icft and a — side lying at 
the right when traversing the contour according to its orientation. 

RH problem 2.1. For e M, we look for a 2 x 2 matrix-valued function v) 
satisfying 

(1) ífíC; v) is analytic for C e C \ S*; 

(2) Í'+(C;í^) = '\!-{C,;v)Jk: for C e Ffc, fc = 1,...,4. Here ^± stands for the 
limiting valué of ^ at the ± side of Vk and Jk is as in Figure [1] 

(3) As C ^ 00 we have 

(4) '^{C,] v) is bounded near C O- 

This RH problem was introduced by Flaschka and Newell in [31]. They 
showed that from this RH problem one can retrieve the Hastings-McLeod solu- 
tion for the Painlevé II equation. More precisely, define q{v) by 

q{v)= hm C*i2(C;i^)e~^^'^'-^'"^, 
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where ^12 is the 12-entry of 'I', thcn q is the unique solution to the Painlevé II 
equation 

q"{y)=2q{vf + yq{y), 

uniquely characterized by the asymptotic condition q{v) = Ai{v){l + o(l)) as 
V — ^ +00. One can show that there exists a unique solution to RH problem 
I2.1l if and only if the Hastings-McLeod solution q has no pole at v. Since it is 
known that this solution has no real polas [37], it follows that exists for all 
€ M. The kernel Kpn at the right-hand side of (|2.8I) is now given by 

(2.9) KMx.y.y) ^ 2ni{l-y) -l) *"'(y; 

where i^) is the unique solution to RH problem 12. II 

For interior singular point of higher ordcr, p{x) ^ {x — a;*)^™ the kernel 
is characterized by a RH problem associated to a special solution to the m-th 
member of the Painlevé II hierarchy. 

2.3.2 Singular edge points 

Near a singular endpoint we have similar limits for the kernel but where the 
Painlevé II equation is replaced by the Painlevé I equation. For example, in [16j 
it is proved that the singular endpoint with vanishing exponent 5/2 is related 
to the second member of the Painlevé I hierarchy. For the general exponent 
m + 1/2 related to the m-th member of the Painlevé I hierarchy see [T4] . 

It is interesting to note that is not possible to obtain a singular edge point 
with an exponent of 3/2 in the Hermitian one matrix model. However, by 
adjusting the model and considering orthogonal polynomials in the complex 
plañe, a vanishing exponent of 3/2 can be reahzed [8j[24l[32] and this singular 
case is related to the Painlevé I equation. However, there is no probabilistic 
interpretation for this situation. 

3 Two matrix model 

In the Hermitian two matrix model we consider the probability measure on the 
space of couples (Mi, M2) of n x n Hermitian matrices given by 

(3.1) ^ exp (-n Tr {V{Mi) + WiA'h) - rAhAh)) áNháh'h. 

Here V and W are two polynomials of even degree and positive leading coefñ- 
cients, Z'^^ is a normalizing constant and r > O is called the coupling constant. 
Note that if t = O, the probability measure factorizes and Mi and M2 are 
independent matrices taken randomly from (|2.ip with potential V and W re- 
spectively. 
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We will be concerned with the asymptotic behavior of the eigenvalues of Mi 
and M2 as n — >■ 00. In this section we will discuss some general characteristics 
of the two matrix model and in particular the relation to certain biorthogonal 
polynomials. 



3.1 Master loop equation 

An important motivation for studying the two matrix model is that there is 
strong evidence that it generates a wide class of singular points that can not 
appear in the one matrix model. 

By formal calculations one can show that the limiting eigenvalue distribu- 
tions are characterized by an algébrale curve, which is also often referred to 
as the master loop equation and is the equivalent of (|2.4|) for the two matrix 
model. To this end, define the following functions 



Yn{x) = V'{x) 



72M 



2M 



Tr 
Tr 



1 



X- Mi_ 
1 

y-M2 



Pn(x,y) = -E: 



2 Ai 



Tr 



V'{x) - r (Mi) W'{v) - W'{M2) 



x- Mi 



y - M2 



En{x, y) = t{V'{x) - Ty){W'{y) - rx) - Pn{x,y) + t\ 

where E^*^ stand for the expectation with respect to p.ip . The conjecture is 
that these functions have expansions in l/n^, i.e. 

y„ = + X;Y^'^ + . . . , 



and similarly for X„,P„ and En. Moreover, the various terms in the expansión 
satisfy a recursive system of equations that are called the loop equations [29] 
(for a survey on the formal analysis of matrix models using loop equations see 
for example [52] and the reference therein). The first of the loop equations, also 
called the master loop equation, reads 



(x, iy(o)(x))=i?(")(ix(o)(y),y)=o, 



(3.2) 



which is to be interpreted as the equivalent of (|2.4p for the two matrix model. 

The function E^^^ is a polynomial of degree deg V in x and degree deg W 
in y, whereas in the one matrix model (j2.4p always has degree 2 in one of the 
variables. Under special cholees of parameters the curve defined by p.2p has 
interesting singular points that cannot occur in the one matrix model [iTj. For 
example, it is possible to obtain a limiting measure where the density vanishes 
with a vanishing exponent p/q. Near these points we expect to obtain new 
interesting scaling limits. We recall that the only possible singular cases for the 
curve (j2.4p corresponding to the one matrix model, are the ones listed in Section 
12.31 (and henee q — 2). 
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At this point we want to emphasize that (I3.2p is derived in a formal way. 
In [35] it was proved that the empirical eigenvalue distributions of Mi and M2 
have weak limit almost surely. This means in particular that Yn and X„ have 
limits, but to the best of my knowledge there is no rigorous proof of an í/n^ 
expansión or of (|3.2|) . 



3.2 Biorthogonal polynomials 

The important feature of the one matrix model that made it possible to analyze 
it explicitly, is that it can be integrated in terms of orthogonal polynomials. In 
the two matrix model there is a similar structure. Let {pj,n)j and {qk,n)k be 
two sequences of monic polynomials, with degpj_„ = j and degqk^n — k, such 
that they satisfy the biorthogonality relation 

P,,„(x)gfc..(y)e-"(^(-)+'^(^)---^)) dxdy = hlj,,, 

for certain constants ^ . Since the orthogonality is not with respect to a 
Hermitian inner product, it is not a priori clear that the polynomials exists. In 
[28] it was proved that they do exist, are unique and have properties that are 
typical for orthogonal polynomials. They have real and simple zeros [28] and 
the zeros satisfy an interlacing property |23| . Their integrable structure has 
been extensively explored in [21 El |3] . 

As in the one matrix model (|2.5I) , the marginal densities or correlation func- 
tion for the eigenvalues have a determinantal structure. To this end, we define 



and the following four kernels 

ri-l 



(3.3) K^^\xi,X2) = ^ -r2—Pk,n(xi)Qk,n{^2), 



— n fe,' 



fc=0 
n-1 



— n fe 



fc=0 

Tl-1 



Kii\y,x) = E T^PkAy)QkAx) - e-"(^(-)+^(^)— 

fe=0 fe." 

n-1 ^ 

-^22^(yi,2/2) = E ■j9—Pk.n{yí)^kAy2)- 



fc=0 fe." 

Then by the Eynard-Mehta theorem [30] the marginal densities of the point 
process given by the eigenvalues of Mi and M2 have the following structure 
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! 2/1 j 



, 2/„)dxfc+i • • • áxnáyi+i • • • dí/„ 



n—k-\-n—l times 



{n-k)\{n-l)\ 



K. 



n'. 



12 



dct 



for k,l = 1, . . . , n. For example, if we average over M2 then we see that the 

(n) 

eigenvalues of Mi form a determinantal point process with kernel K^^ given 
in p.3p . To find the asymptotic behavior of the eigenavalues and the hmiting 
local correlations, it suffices to compute the asymptotic behavior of the kernels 



(") 



3.3 Riemann-Hilbert problem(s) 

The asymptotic behavior of orthogonal polynomials that appear in the one ma- 
trix model can be effectively computed using the RH problem approach (see for 
example 01] for a discussion). It is therefore natural to search for a characteri- 
zation of the biorthogonal polynomials in terms of a RH problem. Several such 
characterizations exist ¡Sj EHl |42l |43] . Here we will discuss the RH problem from 
¡43] as this was the starting point for the analysis in ¡22l [25l ESI |49] . It should 
be noted that this RH problem is equivalent to the RH problem in [5]. 
Write dw — deg W and define the functions 



yjQ-n{W{y)+V(x)-Txy)^y^ 



J =0,l,...d 



w 



Then the kernel k[^^^ in p.3|) can be characterized by the following RH problem. 
RH problem 3.1. We look for a dvv ^ dw matrix valued function Y such that 
1. y is analytic in C \ R: 

/l Wo^n{x) ■■■ Wdw-l,n{x)\ 

O 



2. Y+{x) = Y^{x) 



Vo 



1 / 



for X £ 



as z 00. where rio is 



3. Y{z) = (/ + 0{l/z)) diag(z", z-"", . . . , z""'*"'-^ 
the integer part of {n + dw — j ~ 2)/{dw — !)• 

There exists a unique solution to the RH problem l3.1l which can be expressed 
in terms of the biorthogonal polynomials ¡43] . In particular we have Yii(z) = 
Pn,n{z) and 

('\ 

K''^\xi,X2) = 2^Ux,^x^) (O Wo^n{x2) ■■■ Wrf^ -2 (^2 )) ^ (a;2)> + (a^l ) 
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For the polynomials qn,n and the kernel a similar RH problem holds. 

The strategy for analyzing the two matrix modal, is to perform a steepest 
descent analysis on the RH problem as n — >■ 00. For the general situation this 
is still an open problem. However, for the special case of quartic W we recently 
analyzed the asymptotic behavior in a series of papers that we will discuss in 
the next section. 



4 Two matrix model with a quartic potential 

In this section we will discuss the two matrix model p.l¡) for the following 
special cholee of potentials 

(4.1) V even and Wiy) = i/ + f y^. 

In [23 [5Ü1 |3n] we performed the steepest descent analysis for RH problem 13.11 
for the kernel k[^^ in p.3p that characterizes the eigenvalues of Mi. To this 
end, we used a vector equilibrium problem that we will discuss in Section 4.1 
and 4.2. In Section 4.3-4.5 we further assume that V{x) = ^x'^ and identify 
all possible singular cases that can occur and present a phase diagram [53] in 
Section 4.3. In particular, there is a new type of singular point around we obtain 
a new kernel [55] which we discuss in Section 4.4 and 4.5. 

4.1 Vector equilibrium problem 

Let US assume that we are in the situation í|4.ip . The key ingredient in the 
asymptotic analysis of [55] [5ül |3Ü] is that we found a coulomb gas interpretation 
for the limiting distribution of the eigenvalues of Mi that we will now describe. 
See also [27] for an alternative discussion. 

For two probability measures fj, and v we define the mutual logarithmic 
energy v) and the logarithmic energy by [53] 

= log I ^ I dfi{x)di'{y), = I{fi,fJ.). 

JJ \x-y\ 

The equilibrium problem is to minimize the energy functional 1^2, i's) de- 

fined by 

3 2 

(4.2) i?(i/i, 1/2, 1^3) = Y. ^(^j) - E 1(^3 ^^J+í) 



+ j Viix) diyiix) + J Vsix) di^six) 

among all vectors of measure {vi, 1^2, ^a) satisfying the following conditions: 

1. í^i is a measure on M with total mass 1. 

2. 1/2 is a measure on iM with total mass 2/3. 
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3. 1^3 is a measure on M with total mass 1/3. 

4. V2 < 0'2- 

We need to clarify the external fields Vi and V3 acting on /ii and ^3, and the 
constraint 02 on ^2- To this end, note that the function 

S I— > W(s) — TXS 

with W as in (|4.1[) . has a global minimum attained at some point s ~ s\(x). 
Moreover, if a < O and < |(— a/3)^/^ the function has another local mini- 
mum at some point 52(2;) and a local máximum at some point 53(2;). 
The external field V\ in (|4.2I) is defined as 

Vy{x) = V{x) + W{si{x)) - TXSl{x). 

Moreover, the external field V3 is defined 



W{s3{x)) - TXS3{x) - (W{S2ÍX)) - TXS2{x)), if 52,3(0;) Cxist, 

O, otherwise. 



Finally, the constraint a2 is a measure that is absolutely continuous with respect 
to the Lebesgue measure with density 

'^'i^^^f = — max Res, z £ iR. 

\az\ TT s^+as=TZ 

The following theorem was proved in [2F for the case a — O and in f26^ for 
general a. See also pü], 

Theorem 4.1. Let V and W as in (|4.1¡) . The energy functional E defined in 
(|4.2p has a unique minimizer (/ii, /X2, /X2) among all vectors of measures satisfy- 
ing the conditions Usted below 



The support of the measures in the unique minimizer {fii, (12, fJ-s) have the 
following structure 

^(A^i) = uf=iK-,y 

S{a - ^2) = iM\ (-ÍC2,ÍC2), 

5(aí3) =»\(-C3,C3) 

for some ci,C2 > O and ai < 61 < «2 < . . . < ¿at and N E N. All measures 
are absolutely continuous with analytic densities (possibly except at the origin) . 
Moreover, if C2 > O then the density of 172 — M2 vanishes as a square root at 
±ic2. Similar ly for C3 > O and /Í3. 

Away from the origin, the measure /ii in the minimizer has the same behavior 
that one finds for the equilibrium measure in the one matrix model. Generically 
its density is strictly positive in the interior of the support, it vanishes as a 
square root at the endpoints and the variational inequality is strict. If we are in 
this situation and on top of that we have that in case O € S{a — ^2) or O G S{iiz) 
then the density for that measure is strictly positive at the origin, then we say 
that (V, W,t) is regular. The following result is Theorem 1.4 in ¡26] . 
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Theorem 4.2. Let {V, W, r) be regular. Then, as n ^ oo and n = O mod 3, 

the mean eigenvalue distribution of Mi converges weakly to the first component 
fii of the minimizer (/zi, /Z2, /its) of the vector equilibrium problem. 

We strongly believe that it holds also in all the singular situations and that 
we can drop the condition n = O mod 3 that we needed for technical reasons. 

The measure /ii describes the hmiting behavior for the eigenvalues on the 
global scale and the natural question rises what happens on the local scale. In 
the regular case, the kernel k[^^ converges to the sine kernel in the bulk and 
Airy kernel at the edge points as in (|2.6p and (|2.7p . Away from he origin, the 
measure fii can have the same singular points as in the one matrix model and 

(n) 

the kernel K^^^ has the corresponding limiting behavior. However, due to the 
more complicated interaction with the measure fi2 and /Í3, at the origin new 
critical phenomena may take place which we will explore further in Sections 
4.3-4.5. 

4.2 Associated Riemann surface 

At first sight, the characterization of the limiting eigenvalue distribution in terms 
of a vector equilibrium problem appears to be quite different from the master 
loop equation (j3.2p . However, the solution to the vector equilibrium problem 
can be described by an algébrale curve as we will now show. 
Let TZ = Uj=i be where 

ni = c\s{fii) 

7^2 = C\(5(Aii)U5(a2-//2) 
7^3 - C \ {S{a2 - M2) U S{fi3)) 
7^4 = C \ 5(/i3) 

and TZj is connected to Tíj+i in the usual crosswise manner. The following 
result foUows from the variational conditions for the equilibrium problem and 
is Proposition 4.8 in |26j . 

Proposition 4.3. The function ^ : TZi — >■ C defined by 

(4.3) d/ii(s), 

J z - s 

extends to a meromorphic function on TZ, with a pole at infinity of degree deg V — 
1 071 the sheet TZi and a simple pole at infinity at the other sheets. 

From the last result it follows that the function ^ is described by an algébrale 
curve 

e + P3ix)e + P2{x)e + Pi{x)i + p^{x) = o 

for certain polynomials Pj{x). Moreover, the restriction of ^ to the first sheet 
7^l plays the role of F^o) in 
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Figure 2: The phase diagram in the ar-plane: the critical curves ~ a + 2 
and ar^ = — 1 sepárate the four cases. 

4.3 Phase diagram for quartic/quadratic case 

Let US consider the special situation 

(4.4) Vix) = íx^, and Wiy) ^ + ^y'\ 

and discuss the structure of the vector (/ii , fj,2, fJ-a) minimizing the energy E and 
the associated surface TZ. For this situation an ahernative and perhaps more 
direct derivation of the equiUbrium problem can be found in [23;. 

In this special situation we also have that /ii is supported on one or two 
intervals. Henee there exists a, ci, C2, C3 > O such 

= [-a, a] \ (-ci,ci), 

5(0-2 - ^J■2) = ÍM \ (-ÍC2,ÍC2), 
5(/^3) =K\(-C3,C3). 

If ci =0 the /ii is supported on one interval and if ci > O it is supported on 
two intervals. We now distinguish four different cases 

Case I ci , C3 = O and C2 > O 

Case II C3 = O and ci,C2 > 0. 

Case III C2 = O and ci, C3 > O 

Case IV ci = O and C2, C3 > 0. 

In Figure 13] we showed the sheet structure of the associated Riemann surface TZ 
in each of the four case. Depending on the valúes of r and a we are in one of 
these four cases (or a transition from one to the other). 

In Figure [5] we plotted an ar-phase diagram. The ar-plane is separated into 
four regions by the curves = 2 + a and ar = —1 and each región corresponds 
to a particular case. As long as we are not on one of the two separating curves, 
we are in the generic situation and the local correlation for the eigenvalues of Mi 
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are given by the sine kernel at the bulk and Airy kernel at the edge. However, 
if we pass one of the separating curves there is a transition and at the origin 
the local correlations may have different limits. We will now discuss the various 
transitions but leave the multicritical situation (a, r) = (—1,1) to the next 
section. See also Figure[S]for the Riemann surface TZ at the points of transition. 

Case I ^ Case II. In this transition the support S{iJ,i) splits from one 
interval into two intervals. At the splitting point, the density oí fii vanishes 
quadratically. The situation is the same as in the one matrix model and again 
the local correlations are governed by the kernel pn as given in (|2.9I) . 

Case IVo Case I. Here we see that there is a transition in the support 
of S{fj,3). In Case IV there is a gap in the support 8(^3) that closes at the 
point of transition to dissapear in Case I. In fact, the transition is similar to the 
transition from Case I to Case II, but now it concerns the measure /Í3. As it 
turns out, it does not have an effect on the eigenvalues of Mi. 

Case III ^ Case IV. When we travel from Case III to Case IV the two 
intervals in the support of S{pi) merge to one interval, but simultaneously a gap 
in the support of S{a2 — IJ-2) opens up. The origin is a branch point connecting 
the first three sheets of the Riemann surface. This situation cannot happen in 
the case of the one matrix model. However, it has appeared before in the one 
matrix model with external source [Ü]. The local correlation in this case are 
given by the Pearcey kernel 

(4.5) KM-,y-s) = j^J r ei-^~i»^---i^^+i^-- 

(27ri)^ Je 7-ioo z-w 

where C is a contour that consists of two rays from ±ooe™/^ to O together with 
two rays from O to icx^e"''^/'*. See also pl| [T2l [5T| [57] . 

Case II o Case III. The situation is very similar to the transition from 
Case III to Case IV but now the origin is a branch point connecting the thrce 
sheets 7I2 , TZs and TZ4 . As in the transition from Case IV to Case I it does not 
effect the eigenvalue distribution of Mi. 

4.4 Critical point 

We now deal with the multicritical point (a,r) ~ ( — 1,1). In this case, the 
origin connects all four sheets of the Riemann surface TZ, see also the picture 
at the bottom of Figure [S] In this case the function ^ in (|4.3p is given by the 
algébrale curve 

and the limiting eigenvalue distribution vanishes with a square root near the 
origin which is at the interior of the support. In ¡22j we characterized the 
limiting kernel which we will now discuss. To this end we first need the foUowing 
RH problem. 

RH problem 4.4. Let s,í e R. We search for a 4 x 4 matrix-valued function 
M{() satisfying 
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(1) M is analytic for C G C \ Sa/; 

(2) M+(C) = AÍ_(C)Jfc, for C e Tfc, fc = O, . . . ,9; 

(3) As C — oo with C e C \ Saí we have 

(4.6) Af (C) - (/ + 0{C')) diag ((-C)-'/^ C-'/', C'/') 
/ 1 O -¿ ON 



Vo í o ly 



e3 



2 (_^)3/2+2.(-C)^''^+íC Jc'/'+2<^/^-íC 



(4) Aí(C) is bounded near C = 0. 

The fractfonal powers in C C^^^i C C^^^ ^-nd C ^-^^ C^^^^ are chosen 
such that these maps are analytic in C\ (— oo, 0] and take positive valúes on the 
positive part of the real Une. The fractional powers in C i— ^ (— C)^^^: C ^ (~C)^''^ 
and C ^ (^C)*^^* a-re chosen such that these maps are analytic in C \ [0,oo) 
and take positive valúes on the negative part of the real line. 

The contour T,m is shown in Figure[3]and consists of 10 rays emanating from 
the origin. The function M{í^) makes constant jumps on each of the rays T^. 
These rays are determined by two angles ifi and ip2 satisfying O < (pi < ip2 < 
Tr/2. The half-hnes Tk, fe = O, . . . , 9, are defined by 



and Fs+fe = — Ffe, for fe = O, . . . , 4. All rays are oriented towards infinity. 

For s,í e M the solution to the RH problem exist and is unique [22i Th. 
2.2]. Moreover, the RH problem is related to the Hastings-McLeod solution for 
the Painlevé II equation. Indeed, by rewriting / + 0{í/C)) = I + Aí^^VC + • ■ ■ 
in (|4.6p we have 



i2 



-1/3, 



where the left-hand side is the 14 entry of A/^^^ and q stands for the Hastings- 
McLeod solution. In fact, by taking derivatives with respect to C and the param- 
eters s, t we can obtain from the RH problem a system of first order differential 
equations for M for which the Painlevé II equation appears as the compatibility 
condition pT1[22]. 

We now define Kcr by 



2'KÍ(u — v) 



(1 1 O tí)M{iv]s,tfM{iu-s,ty 



1 
O 

voy 
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Figure 3: The jump contour T,m in the complex C-plane and the constant jump 
matrices Jk on each of the rays Tk, fc = O, . . . , 9. 



where Al"^ stands for the transpose of M and M for the inverse transpose. 
The foUowing theorem is the main result in [35]. 

Theorem 4.5. Let V and W be as in (|4.4p and set 

for a, 6 G M. Then for n — > oo and n = O mod 6, and k[^'^ as in (|3.3p we have 

uniformly for u, v in compact subsets ofM. 

It is interesting that RH problem 14.41 with í = O appeared before [21] in 
the characterization of the kernel near a tacnode singularity in a model of non- 
intersecting brownian paths (see also [Ulll]) for ahernative characterizations) . 
The apparent reason for this is that also in the case of the tacnode singularity, 
the critical measure vanishes as a square root near an interior point. Neverthe- 
less, the kernel describing the tacnode singularity is constructed in a different 
way out of the RH problem. In ¡22) we proved that, perhaps somewhat surpris- 
ingly, the two kernels define essentially different processes. 
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4.5 Reductions of the new kernel 

There is a a certain hierarchy in the limiting kernels that we have discussed so 
far. Let us denote the sine and Airy kernel (i.e. the right-hand sides of (12.61) 
and (|2.7p ) by i^sine and ííTaí- We also recah the definition of Kpu and Kp^ in 
(|2.9p and (|4.5p . The kernel is on top of the hierarchy, in the sense that the 
four other kernels are limit points for Kcr- For instance, from the phase diagram 
in Figure [2] we see that if we walk from the critical point {a,T) = (—1,1) to 
the right along the curve = 2 + a, we end up in the in the critical situation 
described by the kernel Kpu. In j22¡ we proved that Kpu is a limit point of Kcr 
in the following way. There exists a function h such that 

hm 23a-f^iícr 23 ax, 2^ay; ^, -a 1 - ^ = Kp^ix, y; 23a , 
a->-+oo h[y,a) \ 2 \ / J 

uniformly for x, y in compact sets. (Note that as the correlation functions are 
determinantal they are invariant under conjugation of the kernel by the function 
h and the point process is not changed) 

If we walk from {a,T) = (—1, 1) along the curve r = — l/a in the phase 
diagram, then we expect to obtain the Pearcey kernel as a limit. In a recent 
paper [34] the authors proved that 

(4.8) liin —^K„\—f^,—f^,-\a^,a\\-—^\\^Kpc(x,y-,a), 
a^+oo ^2a4 Vv2a4 V2a4 \ 2a2 / / 

uniformly for ce, y in compact subsets. 

Finally, we note that the transitions Kpc Kaí and Kxi — > -K^Sinc can 
for example be proved by classical steepest descent methods on the integral 
representations for the Airy and Pearcey functions. Moreover, the transition 
Kpii — > Ka_í can be found using steepest descent techniques on the RH problem 
[O for Kpii. 



5 Singular points in one matrix models revisited 

Now that we have characterized the correlations for the eigenvalues of Mi it 
is natural to ask about the eigenvalues of M2, especially for the multi-critical 
point in the quartic/quadratic case. We will now answer this question in a 
general fashion and show that the multi-critical point discussed in Section 4 can 
be embedded into a larger class for which a rigorous treatment may be within 
reach. 

First note in case 

V{x) = \x'^ and W arbitrary, 
we can complete the square and rewrite (|3.ip as 

(5.1) ^ exp (-nTr (i(Mi - rMa)^ + WiNh) - ^r^Aía)) dMidMa. 
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Henee the matrices Mi — TM2 and M2 are independent matrices taken randomly 
from the GUE (=Gaussian Unitary Ensemble) and the one matrix modcl (12.11) 
with potential 

Weff(2/) = W{y) - \T^y\ 

Moreover, by writing the matrix Mi as Mi — {Mi — TM2) + TM2, we see that 
Mi is a Unear combination of two random matrices taken independently from 
the GUE and the one matrix model (j2.ip with potential W^ñ respectively. 

It is also important to note that in this situation the spectral curve p.2|) can 
also be computed in a fairly exphcit way, 

(5.2) t(x - tOÍW'^AO + r^e - Tx) ~ QiO = O, 

where Q is the polynomial 

and ¡J^w^íi the equihbrium problem measure minimizing (|2.3p with potential 
Weff ■ We recall that there is no proof that (|5.2p indeed characterizes the limiting 
eigenvalue distribution of the matrix Mi , but we are free to use it as an ansatz 
in a RH analysis and prove this fact a posteriori. This may be a good starting 
point of dealing with some particular interesting examples, some of which we 
will now discuss. 

Let US proceed by choosing the parameters in the potential W so that WcS 
does not depend on r and consider the family of random matrices 

M{t) = Mgue + tMw^,¡ 

parametrized by r. Note that {j^M{t))t interpolates between Mque for t = O 
and Mw^a for t — > 00. Let us also assume that Wcff leads to an equilibrium 
measure /íw^h for which the density vanishes with exponent 2m at the origin. 
Then on the local scale we expect to see a transition from the sine kernel to 
the kernel associated to the singular interior point (related to the Painlevé II 
hierarchy) . 

To see how this transition comes about we rewrite (|5.2p to 

as ^ — > O, for some positive constant c > 0. Set — -y/— W^'^(0)/2. If r 7^ r^r 
then for a; = O we find that ^ = O is a double solution and we have a singular 
point that connects two sheets of the Riemann surface. If the point is on the 
physical sheet TZi then the local correlations are given by the same kernel as 
we have for M\\r^^^. In case r = Tcr then for x = Q the solution ^ = O is 
of multiplicity at least four (in case m = 1) or at least six (in case m > 2) 
and henee the singular point connects several sheets. Here we expect a more 
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complicated kernel, characterized by a RH problem that has the same size as 
the number of sheets involved. 

Concluding, by adding an independent GUE matrix to a matrix taken ran- 
domly from a one matrix model with a singular interior point, we can construct 
higher order critical phenomena. It is an interesting open problem to analyze 
these cases explicitly and identify the local correlations (and their reductions). 



Example 1. In the first example we return to the case in Scction r4.4l and consider 

(5.3) V{x)^\x\ and W {y) ^ - ^-=^y^ 

so we are are on the curve a = — 2 + in the phase diagram in Figure [2] Henee 
we have Wcñ{y) = y'^/^ — y^ and henee, in a double scaling limit, the local 
correlation for the eigenvalues for M2 are governed by the Painlevé II kernel 
(|2.9p . As for the eigenvalues of Mi , from the phase diagram we read off that for 
T < 1 the local correlation are given by the sine kernel and for r > 1 we obtain 
the Painlevé II kernel. The transition takes place at r = 1 (and henee a — — 1) 
where the kernel has the scaling limit (|4.7p . 

Example 2. In the second example we consider the 

W{y)^\y'-^y' + \{-^ + T^)y^ 
For this potential the equilibrium measure for Woff is given by 

ájJL 



da; ^ 

and we have a quartic vanishing at the origin. In this case, the spectral curve 
(15.31) takes the form 



2-4/3^2 ^ ^4 ^ _ ^ ^1/3^ _ ^2^ ^ 2^/3^3 - e) = O, 

In particular, for r = 2^^/^ -^q have 

x"^ + f - xf{2^/^^^ - 2) = O, 

and from here we see that for a; = O the solution ^ = O has multiplicity six and 
the origin connects all sheets. 

Example 3. In the final example we show that similar phenomena occur when 
dealing with singular edge points. Consider the potential 

W{y)^ly + {\ + ^4)y^-±y + ^,y\ 

In this case, the effective potential W^ñ gives rise to an equilibrium measure 
with a singular endpoint 

^ = ^(.-2)^v/ÍT^. 
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See also [T^. For thc special case t — \fhl^ the spectral curve takes the form 

(-32 + 16y - 8?/2 + 4í/3 _ y^) + ^5(8 + 4y - 4y2 + y^)x - 5a;^ = O 

For a; = 4-\/5/5 we find four solutions ^ = 2 and henee this point connects all 
sheets. 
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Figure 4: The Riemann surface in the four diffcrcnt cases that are possible in 
the case V{x) = ^x'^ and W{y) = \y'^ + ^y^. 
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Figure 5: The Riemann surface in the transition from one case to the other for 
V[x) = and W{y) = + f y^. 
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